Quantum fidelity in the thermodynamic limit 
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We study quantum fidelity, the overlap between two ground states of a many-body system, fo- 
cusing on the thermodynamic regime. We show how drop of fidelity near a critical point encodes 
universal information about a quantum phase transition. Our general scaling results are illustrated 
in the quantum Ising chain for which a remarkably simple expression for fidelity is found. 
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A quantum phase transition (QPT) happens when dra- 
matic changes in the ground state properties of a quan- 
tum system can be induced by a tiny variation of an 
external parameter This external parameter can be 
a strength of a magnetic field in spin systems (e.g. Ising 
chains [2[ and spin-1 Bose-Einstein condensates [3|), in- 
tensity of a laser beam creating a lattice for cold atom 
emulators of Hubbard models 4(, or dopant concentra- 
tion in high-Tc superconductors [5|. At the heart of the 
sharp transition lies non-analyticity of the ground state 
wave- function across the critical point separating the two 
phases. QPTs, traditionally associated with condensed 
matter physics, are nowadays intensively studied from 
the quantum information perspective (see e.g. Q). 

Quantum fidelity - also referred to as fidelity - is a 
popular concept of quantum information science defined 
here as the overlap between two quantum states 



F(g,6) = \(g-6\g + 6)\, 



(1) 



where \g) is a ground state wave- function of a many-body 
Hamiltonian H(g) describing the system exposed to an 
external field <?, and 5 is a small parameter difference. It 
provides the most basic probe into the dramatic change 
of the wave- function near and at the critical point Q • 

The recent surge in studies of fidelity follows discov- 
ery that quantum criticality promotes decay of fidelity 
Q. This is in agreement with the intuitive picture of 
a QPT: as system properties change dramatically in the 
neighborhood of the critical point, ground state wave- 
function taken at different values of the external param- 
eter - \g — 5} and \g + 5) - have little in common and so 
their overlap decreases. 

As fidelity is given by the angle between two vectors in 
the Hilbert space, it is a geometric quantity Q. Thus, it 
has been proposed as a robust geometric probe of quan- 
tum criticality applicable to all systems undergoing a 
QPT regardless of their symmetries and order parame- 
ters whose prior knowledge is required in traditional ap- 
proaches to QPTs. Fidelity has been recently studied 
in this context in several models of condensed matter 
physics (see [§] and references therein). 

Besides being an efficient probe of quantum critical- 
ity, fidelity appears in numerous problems in quantum 



physics. Indeed, it is related to density of topological 
defects after a quench 1^- 13] , decoherence rate of a test 
qubit interacting with an out-of-equilibrium environment 
[l3| , orthogonality catastrophe of condensed matter sys- 
tems (see 14| and the references citing it). Therefore its 
understanding has an interdisciplinary impact. 

To unravel the influence of quantum criticality on fi- 
delity one has to determine if its drop near the critical 
point encodes universal information about the transition 
in addition to providing the location of the critical point. 
This universal information is given by the critical expo- 
nents and reflects symmetries of the model rather than its 
microscopic details. In the "small system limit", which 
broadly speaking corresponds to 5 —> at fixed system 
size N, the answer is positive. This is ex plor ed in the 
fidelity susceptibility approach where @, 0, KB] 



F(g,6)~l-6 2 XF (g)/2, 



(2) 



and xf stands for fidelity susceptibility. Universal infor- 
mation, or simply the critical exponent v, is encoded in 
its scaling: at the critical point Xf(#c) ~ N 2 / dv , while 



N\g-g c 



where d is system 



far away from it yp 
dimensionality [ill Qjl lI6( . 

In the thermodynamic limit, which broadly speaking 
corresponds to N — > oo at fixed S, the answer is positive 
as well. This is our key result stating that 



(3) 



where A is a scaling function. In particular, we see 
that fidelity is non- analytic in 5 at the critical point, 
\nF{g c ,5) ~ — N\8\ dw , while away from it, i.e., for 
\8\ <C |.g — g c \ *C 1, we obtain 



lnT(g,S) ~-N6 2 \g- gi 



\dv-2 



(4) 



after expansion of the scaling function. These results, 
in particular, set firm foundations for usage of fidelity 
as a probe of quantum criticality in thermodynamically- 
large systems. In the context of theoretical studies of 
QPTs, the strength of the fidelity approach lies in its 
simplicity: all information encoded in the ground state 
wave-function(s) is "compactified" into a single number. 
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FIG. 1: (color online) Fidelity of the Ising chain near the 
critical point as a function of (a) the system size TV at fixed 
S = 1CP 4 and (b) parameter difference S at the fixed system 
size TV — 10° . On both panels the curves from top to bottom 
correspond to J-(l, 5), J-(l + 5, S) and T(l + 5<5, 5). 



A competing approach for extraction of the exponent v 
- study of the asymptotic decay of correlation functions 
to obtain the correlation length - is considerably more 
complicated. Below we illustrate these predictions on 
the paradigmatic model of quantum phase transitions, 
the Ising chain, and discuss the scaling theory that leads 
to © and (J). 

The Hamiltonian of the one dimensional Ising chain 
reads {h] 



H(g) 



N 

■E< 

i=l 



where g stands for a magnetic field acting along the z 
direction. Above the spin-spin interactions try to enforce 
±x polarization of spins, while the magnetic field tries to 
polarize spins along its direction (+z for g > 0). This 
competition results in two critical points at g c = ±1: the 
system is in the ferromagnetic (paramagnetic) phase for 
— 1 < g < 1 (|g| > 1). The critical exponent v = 1. This 
model is solved by mapping spins onto non-interacting 
fermions via the Jordan- Wigner transformation [l[ . 

Behavior of fidelity ([T]) around the critical point, g « 
<? c , is summarized in Fig. [T] In Fig. [TJl the parameter 
difference S is kept fixed and the system size is increased. 
For small system sizes we reproduce the known result, 
In J 7 ~ —7V 2 [ 7L re sulting from finite size scaling effects 
(see e.g. (9HT1L 12. Tfjj]). For large system sizes, however, 
we obtain In J 7 ~ —TV in agree ment with © and the 
fidelity per site approach 17Hl9{. As is shown in Fig. [2] 
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FIG. 2: (color online) Study of the crossover between the 
"small system limit" and the thermodynamic limit illustrated 
in Fig. [T] As the system size is increased in Fig. the 
slope of the curves changes smoothly from 2 (corresponding 
to lnJ-' ~ —TV 2 ) to 1 (corresponding to InT ~ —TV). The 
crossover region between the two limits is centered around 
TV = TV3/2 where the slope equals 3/2. To find it, we have 
calculated numerically J-(g, S) vs. TV - as in Fig. [jj, - for 
various 5's and found that TV 3 / 2 |<5| ~ 1- This is illustrated 
in this figure where data sets from top to bottom correspond 
to results obtained for g — 1, 1 + S and 1 + 5<5, respectively 
(similarly as in Fig. []J,). The power-law fits (straight lines) 
to numerical data (crosses) give TV 3 / 2 = a 1 5 1 ~ b , where b — 
0.995 ± 0.003 for all three fits, while the prefactor a changes 
between the fits from 3.6 to 0.3. Similar analysis can be done 
on curves shown in Fig. [lb providing the same result. Thus 
we conclude that the crossover condition reads TV|<5| ~ 1 near 
the Ising critical point. 



which will be discussed below. 

Similarly, we observe two distinct regimes when the 
system size TV is kept fixed and the parameter differ- 
ence S is varied (Fig. [TJ)). For TV|<5| « 1 we observe 
In J 7 ~ — 5 2 , in agreement with @, while for TV|i5| 1 
we find In J 7 ~ — \S\ in agreement with ([3]). In the lat- 
ter fidelity approaches non-analytic limit (where dgJ- at 
S = is undefined) reflecting singularities of the ground 
state wave-function resulting from the QPT [2(| • 

We also see on both panels of Fig. [1] that all curves 
collapse for N\8\ <C 1, while they stay distinct in the op- 
posite limit. Thus, for TV|i5| 3> 1 sensitivity of fidelity to 
quantum criticality is enhanced. This can be understood 
if we focus on Fig. QJi: in the large TV limit dramatic 
changes in the ground state wave-function near the crit- 
ical point are expected. 

As analytical results for fidelity are scarce, we find it 
remarkable that we can derive accurate analytical de- 
scription in the complicated limit of N\5\ ^> 1, where 
the Taylor expansion ([2]) fails. To proceed, we calcu- 
late J-"(l + e,<5), where e measures distance from the 
critical point. For the Ising chain J- = ^ik>ofk, 
where f k = cos(6»+(fc)/2 - 6L(fc)/2) and tan(0±(fc)) = 
sin£;/(l + e ± S — cos A;). We stay close to the criti- 
cal point so that < |e| <C 1 and introduce nat- 
ural parameterization: c = e/|<5|. Taking the limit of 
TV — ^ 00 at fixed S the product Ilfc/fc can be changed into 
exp(TV J dk In fk/2n), which can be further simplified to 



the transition between the two regimes takes place when 



In J" 



N\S\A(c) 



(G) 



N\5\~l, 



(5) 



in the leading order in S and e. This result is in prefect 
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FIG. 3: (color online) Fidelity J-(g, 5) of the Ising chain 
near the critical point: thermodynamic limit (main plot) vs. 
"small system limit" (inset). Main plot: black curve is our 
analytic approximation ([6]), while red crosses come from nu- 
merics. Both were obtained for N = 2 x 10 5 and 5 = 10~ 4 
(N\S\ > 1). Inset: numerical result for N = 10 3 and 5 = 10~ 4 
(N\8\ <C 1). In the "small system limit" fidelity stays close to 
unity at any distance from the critical point, while in the ther- 
modynamic limit it can interpolate between zero and unity. 



agreement with our universal scaling law @: note that 
v, d = 1 in our model and c = (g — <? c )/|<5|. Moreover, 
it agrees well with exact numerical simulations: Fig. [3] 
Above A(c) is given by 



A(c) 
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|c| < 1 



|c| > 1. 



(7) 

where c x = -4|c|/(|c| - l) 2 , c 2 = (\c\ + l) 2 /(|c| - l) 2 , and 
K and E are complete elliptic integrals of the first and 
second kind, respectively. Agreement between ([7]) and 
numerics is very good: see Fig. |4]for detailed comparison 
of A{c) to numerics. Several interesting results can be 
obtained now. 

First, Eq. ([6]) shows analytically how the so-called An- 
derson catastrophe - disappearance of the overlap be- 
tween distinct ground states of an infinitely large many- 
body quantum system [3] ~ happens in the Ising chain. 

Second, Eq. ([6]) explains the lack of collapse of the var- 
ious curves providing fidelity around the critical point in 
the N\S\ ^> 1 limit. Indeed, fidelity calculated for two 
ground states symmetrically around the critical point 
is ~ exp(— N\S\/4:), but if one of the ground 

states is obtained at the critical point, F(\ ±6,5) ~ 
exp(-A r |^|(7r-2)/47r). In the opposite limit of N\S\ < 1, 
T ~ 1 — 5 2 N 2 / 16 in both cases explaining the collapse 
of all curves in this limit in Fig. [TJ 

Third, there is a singularity in the derivative of fidelity 
when one of the states is calculated at the critical point: 
dT(g ± 5,5) / dg\ g= g a= \ is divergent when N — > oo at 
fixed 5. This reflects singularity of the wave-function 
at the critical point approached in the thermodynamic 
limit. Quantitatively, dA(c) / dc\ c ^i± — In 1 1 — c|/4-7r — 
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FIG. 4: (color online) Upper plot: scaling function A(c) of 
the Ising chain. The solid black line provides the analytic re- 
sult (O, while the red crosses show numerics (i.e., In ^?-"/7V | <5 1 ) . 
The inset highlights singularity at c = 1. Lower plot: loga- 
rithmic divergence of dAjdc\ c= \ discussed in the text. The 
solid black line is the derivative of 0, while the red dashed 
line is a numerical result: the difference between the two near 
the pinch point at c = 1 is due to the finite system size N [19l ] . 
It disappears in the limit of N — ► oo. In both plots numerics 
is done for TV = 10 5 and 5 = ttIO" 3 



31n2/47r+(l±l)/8 + 0((l-c) In |1 — c|), which is loga- 
rithmically divergent at c = 1 (Fig. [2}. This divergence 
is a signature of a pinch point found in [l7l - [l9j when fi- 
delity between two distinct ground states states was stud- 
ied. The logarithmic di verg ence in the Ising chain was 
numerically observed in [19j. 

Last but not least, we obtain from © a compact ex- 
pression for fidelity away from the critical point. Taking 
|c| > 1 (but still |e| = \c5\ < 1), A(c) ~ l/16|c| and so 



T 



exp(-iV(5 2 /16|e|), 



(8) 



in agreement with (|4]) . This reduces to a known result for 
fidelity susceptibility when the argument of the exponent 
is small and so T ~ 1 — <5 2 iV/16|e| (see e.g. [9|), but 
provides a new result in the opposite limit where lowest 
order of the Taylor expansion is insufficient. We notice 
also that (jSJ is analytical in 5 even in the limit of N —> oo: 
there are no singularities expected when the system is far 
away from the critical point. 

Below we derive general scaling results ([3]) and (|4|). 
This can be done by studying the scaling parameter 

d(g + 5,g-5) = - lim In F(g, 5) /N, 

introduced in [17] in the context of fidelity per site ap- 
proach to the thermodynamic limit. We expect that this 
limit is reached when 



min[£( 5 + 6), £{g-8)] « L, 



(9) 



4 



where is the correlation length at magnetic field g 
and L is the linear size of the system (N = L d for a 
rf-dimensional system). Indeed, the smaller of the two 
correlation lengths sets the scale on which the states en- 
tering fidelity "monitor" each other ([1]). In particular, 
it explains our results showing that the thermodynamic 
limit is reached even when one of the states is calculated 
at the critical point and so its correlation length is infi- 
nite. Near a critical point (|9]) is equivalent to L\S\" » 1 
[2lj ]. For the Ising chain studied above it reads N\S\ J> 1 
properly predicting the crossover condition §5§ obtained 
from numerical simulations (Fig. [2]). 

Generalizing the scaling theory of second order QPTs 
(Sec. 1.4 of [22]); we propose the following scaling ansatz 
for the universal part of the scaling parameter 

d(g c + e + S,g c + e-S) = b- d f((e + 8)b x ' v , (e - S)^"), 

where / is the scaling function, b is the scaling factor, 
and v is the critical exponent providing divergence of the 
coherence length £ ~ \g — g c \~ v ■ The scaling function 
depends on both e + 8 and e — 8 as they are renormalized 
simultaneously. The factor b~ d appears for dimensional 
reasons. Scaling of e + 8 and e — 8 is given by scaling of 
the correlation length f (e ±S) = 6£((e ± S)b 1 / 1 '). 

Taking g = g c + e, introducing natural parameteri- 
zation e = c\S\, and fixing the scale of renormalization 
through {Sib 1 /" = 1 we obtain d(g + S,g-S) = \5\ du f(c + 
1, c — 1). It gives after setting /(c + 1, c — 1) = A{c). 
In a general context, ([3]) shows how universal part of the 
scaling parameter causes the Anderson catastrophe near 
a critical point. 

The scaling function A(c) can be simplified away from 
the critical point. We assume below e, 8 > for simplic- 
ity, take 5 « e « 1, and set 6 through (e + 8)b 1/u = 
1 exploring the freedom to choose the renormalization 
scale. Simple calculation results in d(g + 8,g — 8) = 
(e + 8) dlJ f{l, (e — 8)/(e + 8)), where the second argu- 
ment of / is close to unity. Expanding / in it we get 
d(g + S,g-S) « 28 2 e d »~ 2 f"{l,x)\ x=1 as f(l,x) has a 
minimum equal to zero at i = 1. Thus, away from a 
critical point we end up with (0}. When the system is 



small enough, NS 



<C 1, but still in the thermody- 



namic limit ([§]) , we reproduce the known result for fidelity 
susceptibility 1 - J 7 ~ 8 2 N\e\ du ~ 2 [ill, 0, El • Otherwise, 
(HJ provides a new result. 

On general grounds, one can expect that for systems 
with dv > 2 non-universal (system-specific) corrections 
to the above scaling relations may be significant [l2j ]. 
which requires further investigation. 

Summarizing, our work characterizes fidelity - a mod- 
ern probe of quantum criticality - in the thermodynamic 
limit. We have derived, and verified on a specific model, 



new universal scaling properties of fidelity. These find- 
ings should be experimentally relevant as the first exper- 
imental studies of ground state fidelity have been already 
done H3. 

This work is supported by U.S. Department of Energy 
through the LANL/LDRD Program. 
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